Cosmological Equation of State and Interacting 

Energies. 

Alberto C. Balfagon(l) and Raul Ramrrez-Satorras(2). 
Section de Ffsica Teorica y Aplicada. 
Departamento de Ingenieria Industrial, 
Instituto Qumiico de Sarria, Barcelona (Spain) 
Via Augusta 390, Barcelona, Spain 

(1) albert.balfagon@iqs.edu 
(2) raulramirezs@iqs.es 



Abstract 

In this paper we study a model of cosmic evolution, assuming that 
the different components of the universe could interact between them any 
time. An effective equation of state (EOS) for the universe has been used 
as well. A particular function for w, which gives a good agreement be- 
tween our results and the experimental data, has been studied. Finally, 
the model obtained has been applied to different important cases. 
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1 Introduction 



At the end of the last century, it was experimentally showed that our universe 
has recently started an acceleration era [1, 2, 3]. One of the explanations to 
this phenomenon is based on the existence of a dark energy, which draw new 
attention to the cosmological constant, see [4,5] for a review and [6] for the 
study of the properties of dark energy from recent observations. This simple 
explanation had some problems that were faced using different models [7, 8, 9]. 

Most of the models assume that the components of the universe do not in- 
teract and look for an equation of state for the dark energy that could solve 
the problems [10], several authors imposed cosmological [11] or energetic [12] 
conditions in order to delimit the parameters. Sometimes, the studied models 
allow the variation of the equation of state with time [13]; there are also models 
that consider the possibility of the evolution of the dark energy density along 
the time [14, 15, 16]. On the other hand, other authors directly change the 
Fricdmann's equations [17]. Several authors investigated the possibility of in- 
teraction, but limited to some components of the universe, for instance [18, 19, 
20, 21] consider the interaction between dark energy and dark matter. Recently 
some authors studied observational constraints for interactions between dark 
energy and dark matter [22,23,24,25]. Phantom cosmological models have also 
been proposed but, aside its physically odd negative energy, they seem not help 
to explain the problem correctly [26,27]. There are very few models that study 
the possibility of an effective EOS for the universe [13, 20], however they only 
work with some components and with red shift not very far from the present. 

In this paper a general model for the universe has been studied, with very 
few restrictions. Thus the number of the constituents of the universe has been 
not fixed and interactions between any kinds of components in any time of the 
history of the universe have been permitted. The work has been carried out with 
an effective EOS for the universe that is taken as p — wp, where any variable can 
depend on z. The paper is divided in three sections; in the first one, the general 
equations that have been used are presented. In the next section a particular 
form of w has been studied, fixing the different parameters in order to get the 
best agreement with the experimental data [1]. Finally, in the last section, the 
model obtained has been applied to some particular situations that let to some 
interesting deductions about the past and future cosmic evolution. 

2 General equations. 

The units used throughout this paper will be physical units in order to be more 
pedagogical and obtain equations more intuitive. As is currently done nowadays, 
the study has been carried out considering that our universe is flat. Friedmann's 
equations with the last assumptions are [28]: 




(1) 
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Where p' stands for any kind of energy but the one that comes from A (the 
so-called cosmological constant, however in this paper, the possibility that it 
could change with z is been taken into account). 

Where p', as in Eq.(l), takes into account any pressure except for the one 
due to A. This pressure p' can be thought as the pressure made by all the 
components of the universe (without the pressure from A ) together with the 
pressure due to all kinds of interactions between all the constituents of the 
universe (considering now the interactions where A appears). Eq.(l) and (2) 
can now be rewritten as: 

Where p is the total energy density of the universe. 

Where p is the total pressure of the universe. From now on, all Friedmann's 
equations that will be used are Eq.(3) and (4), with the last interpretations for 
p and p. 

From Eq.(3) and (4) we can get: 

P = -3^(p + p) (5) 
The cosmic acceleration can be deduced from -: 

a 

- = -(3 P + p)^ (6) 
The relation between the energy density and the pressure will be taken as: 



P = wp (7) 
Where w and p can be functions of the red shift. Applying Eq.(7) in (6) it 



can be deduced the necessary condition for the acceleration of the universe 

1 

3 

Using (5) and (7) we get: 



w < -\ (8) 



p = p Q exp ^3 J Y^ d ^j (9) 

Where sub index stands for the present time. The Hubble function, de- 
duced from (3) and (9), is: 
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The fitting of the different parameters used will be done through the ex- 
perimental data from the SNela [1]. The experimental data are expressed by 
the distance modulus for different red shifts, and its relation with the Hubble 
function is [7,28]: 

D L = (l + z)cJ o - (11) 

Where Dl is the light distance. The distance modulus A is defined as the 
difference between the apparent and absolute magnitudes for each one of the 
SNela supernovas, its relation with the light distance is: 

A = 25 + 51ogL> L (12) 
Where Dl is expressed in Mpc. 



3 Equation of state (EOS) 

In this section one of the many possible forms for the equation of state (7) will 
be studied. The reason for the selection done is similar to the followed by other 
authors [13]: w must fit the EOS for radiation for the early universe, vary with 
time (or as in an equivalent manner, vary with the red shift) and allow that in 
the present time, the universe could attain an acceleration era. The w will be: 

1 (! + •*) (P- 1 ) 

w = b^T7)-^WT7) (13) 

Where a, (3 and 7 are three constants that will be fixed using the experi- 
mental data. With Eq.(13) the Eq.(9), (10) and (12) can be expressed as: 

z + (3\ 3 '' (a + z) M t x3(i- 7 ) 



<-"\rr) ^<i+.r-» (») 



P J V a 



c(l + z)(a[3^) 1/2 f z 1 







A = 25 + 51og^ ^ ' I — — s<i=wdz (16) 



{z + 13)- (a + z)* (l + zY 



The limits studied for each one of the parameters a, (3 and 7 have been: 
ae [3,4.5],/? G [3, 4.5], 7 G [0.5,1.5] 
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These limits were fixed in order to let w fit, as much as possible, the known 
history of the universe, mainly the recently accelerated era for z between 1 
and 2. To study the best set of parameters a, (3 and 7 a full factorial design 
of experiments (DOE) [29] has been implemented, in order to obtain the best 
fitting between Eq.(16) and the 307 SNela distance modulus experimental data 
[1] Taking Ha = lOO/i , where h can be thought as a new free parameter to be 
fixed, the best set of parameters were: 



a = 4.1,(3 = 3.5,7 = 1-0,^= 0.697 



(17) 



With the set of parameters above mentioned (17) the mean absolute error 
found was 0.23, while the mean absolute error for a model that consider no 
interactions and only with baryonic matter (w = 0), dark matter (w = 0) and 
dark energy (w = —1), was 0.23. Eq.(13)-(16) can now be expressed, taking 
7 = 1.0, as: 



w 



l (l + z) Q3-1) 
3 {a + z) {(3 + z) 



P = Po 



z + /3\ (a + z) 



H = Hn 



z + /3 
(3 



3/2 



a + z 



1/2 



A = 25 + 5 log 



,c(l + z) (a/? 3 ) 
Ho (a - 13) 



1/2 



(3 + z 



(18) 
(19) 
(20) 
(21) 



In Fig.l it can be seen the evolution of Eq.(18) where z G [—1,10]. The 
red shift for the start of the acceleration (8) is inside z € [1.6, 1.7]. It can be 
deduced that the EOS of the universe will tend to a cosmological constant EOS: 
w — > —1, then p — > —p. 



w(z) 




Figure 1: Function w(z) against z, using constants(17) 
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In Fig. 2 the evolution of the total pressure (7) in function of z, within the 
interval z E [—1,8] it is shown. 
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Figure 2: Function p(z) against z, using constants (17) 

In Fig. 2 it can be observed that the universe reached a total negative pressure 
era, before the present acceleration era, where gravitation still dominated during 
an interval of time. The continued expansion of the universe made possible that 
the negative pressure overcame gravitation driving to the accelerated cosmic 
expansion. 

In Fig. 3 the experimental values [1] are plotted against the obtained from 
Eq.(21) using the parameters in (17). 



A 




0.5 1.0 1.5 



Figure 3: Distance modulus against z. Points are experimental data and the 
line represent the values calculated with (21) using (17). 

Figs. 4, 5 and 6 are magnifications of the graphs showed in Fig. 3 for z E [0, 0.5] 
Fig.4, z E [0.5, 1] Fig.5 and z E [1, 1.5] Fig. 6. 
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Figure 4: Magnification of Fig. 3 for z between and 0.5. 
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Figure 5: Magnification of Fig. 3 for z between 0.5 and 1 
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Figure 6: Magnification of Fig. 3 for z between 1 and 1.5. 



An important issue is to know, as much as possible, the evolution of the 
cosmic acceleration. Using Eq.(6), and with some sensible assumptions, it is 
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possible to get some useful information about the evolution of the cosmic accel- 
eration. Assuming that the scale factor a is an ever growing function (reason- 
able assumption with the present experimental data that we have) , if we express 
Eq.(6), using (7), as: 

then if — (3w + l)p raises this will imply that the cosmic acceleration will 
raise, if it diminishes we cannot obtain any conclusive conclusion because what 
could happen is that the acceleration diminishes or raises more slowly than the 
scale factor a. In Fig. 7 we have plotted — (3w + \)p in function of z. 



-(3w+l)p 
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Figure 7: Representation of — (3w + l)p in function of z. 



From the data (as it can be seen in Fig. 7) it can be deduced that the cosmic 
acceleration rose until z 6 [0.31,0.33] and from this point forward it is not 
possible to assure the exact acceleration's evolution. Of course, it can also be 
seen that at the point when the cosmic acceleration is positive it will remain 
positive forever. 



4 Particular cases. 

In this section, three models have been studied applying the equations that have 
been deduced in the previous section. 

4.1 Model with interaction between the components of 
the universe and constant values for the dark energy 
density and pressure. 

In this model Eq.(7) is taken as the sum of the next terms: 

p = wp = p m + p r + PA + Pi (22) 
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Where p m is the baryonic and dark matter pressure, p r is the radiation 
pressure, p\ is the dark energy pressure and pi is a pressure that comes out 
from the interaction of the components of the universe. The EOS for each one 
of the components will be taken as if there was no interaction between them, so 
loading any kind of interaction in pi . The pressures for the different components 
are: 

Baryonic and dark matter pressure: p m = 
Dark energy pressure: 



Pa = ~PA = 



c 4 A 



Radiation pressure: 



f 



f 



Pr 



Pr = ^Pro (1 + z) 



(23) 



(24) 



3' 3 r 

From the currently accepted adimensional densities values [28] fl A = 0.7 and 



= 0.0001 the expressions (23) and (24) result on: 

7 3c 2 H 2 



PA 



-0.7- 



IttG 



Pr 



0.0001 3c 2 H§ ,„ , 4 

° 0- + z) 



(25) 



(26) 



3 8ttG 

The interaction pressure can be deduced substituting (18), (19), (25) and 
(26) in (22): 



3c 2 iJ 2 



8nG V P 
3c 2 ff 2 
8nG 



+ py (a + z) (\ (l + z) {(3-1) 



0.7 



a \3 (a + z) 
0.0001 



03- 



+ 



(27) 



In Fig. 8 the interaction pressure for z 6 [— 1, 8] is shown. 



P [Pa] 




Figure 8: Interaction pressure as a function of z 
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From the figure above, different interpretations can be drawn; one of these, 
could be that the interactions between the different components of the universe 
gave a positive pressure while the baryonic matter and radiation densities were 
large. With the expansion of the universe and the decrease of densities, an 
interaction between dark energy and the baryonic matter started to dominate 
(so dark energy would not be as dark as it is assumed to be), which resulted in 
a negative interaction pressure. Eventually, when the baryonic matter density 
was very small, an initially small interaction between dark energy and dark 
matter started to dominate, which resulted in a positive interaction pressure. 

Thus, some experimental sign of an interaction between dark energy and 
baryonic matter could be a signal that the above interpretation could be possi- 
ble. 

Other possible interpretation could be the possibility of a self-interacting 
dark matter [30] that could provide an accelerated expansion phase as it was 
shown in [31, 32] 

4.2 Model without interaction and with a variable dark 
energy pressure. 

In this model Eq.(22) can be put as: 

P = Wp = p m +p r +p A (28) 

The dark energy pressure can be expressed in a general form as: 

PA = wapa (29) 

Working in a similar way as it was done to get Eq(27) the following result 
is obtained: 

3c 2 H2 fz + f3\ 3 (a + z\ (\ (1 + 2) /3-1 

PA= ' 11 ' ! 



8ttG V J V a J\3(a + z) (3 + Z/ 

oV0.0001 (1 A (3o) 



3c 2 ( 0.0001 
8irG 



At this point, different ways can be followed depending on which term in 
(29) it is considered not constant. 

4.2.1 p\ = cte 

From equations (29) and (30), and taking Q A = 0.7, it can be deduced that wa 
will be expressed as: 



1 

WA 



0.7 



z + [3\ 3 fa + z\ fl (1 + z) /3-l\ 0.0001 , , 4 

' — 5 — (! + 2) 



(3 J V a ) \3{a + z) [3 + z 



(31) 
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In Fig. 9 w\ has been plotted, Eq.(31), for z e [—1,8] 




Figure 9: wa as a function of z. 



Comparing Fig. 9 with Fig.l it can be observed that while the w of the 
universe tends to -1, the dark energy wa tend to -0.4 and attain values lower 
than -1 during some time. This difference in behavior appears due to the fact 
that in Fig.l w and the energy density of the universe are free to be time (or z) 
functions, while in Fig. 9 a constant dark energy density is fixed. This result 
is interesting because it shows that when some authors claim for the necessity 
of a phantom period for the EOS of the dark energy, in fact it is not necessary 
because what happen are a variation of the energy density and the w of the 
universe, keeping always w > —1. In the same way, the results obtained show 
that it is very difficult to consider that the dark energy density is constant (it 
leads to a phantom energy), so it must change during the cosmic evolution. 

4.2.2 p A ^ cte 

In this case the EOS for the dark energy density will be: 

Pa = -Pa (32) 

With (32) the Eq.(28) let the determination of the density parameter J7a 
(which nowadays is considered to be Sl A — 0.7) as a function of z: 

_ 0.0001 (1 + z) 4 f z + \ 3 f a + z \ /l (1 + z) 0-1 
A ~ 3 \ ) \ a ) \3(a + z) + z 

In Fig. 10 the Eq.(33) for z € [—1,10] is represented. It can be observed 
that the density parameter took negative values in the far past and nowadays 
(z = 0) tends to 0.63. 
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Figure 10: Parameter density for z running between -1 and 10. 

A negative energy is not very acceptable, so model 4.2.2 is not physically 
correct. However, given the fact that for z = 0a value similar to the accepted 
one is obtained, it can be deduced that if the dark energy density varies with z 
then this variation must be similar to the variation showed in Eq.(33). In order 
to avoid negative values for the density it has to be assumed that w\ can vary 
and/or that there must have an interaction pressure. 

4.3 Model with interaction and a variable dark energy 
pressure. 

In this last case a model with interaction between the different components of 
the universe (expressed as an interaction pressure) is studied. We will also admit 
that the dark energy density could vary, assuming an EOS for dark energy of 
the cosmological constant type: p\ — — pa- The variation of the dark energy 
density will be deduced from Eq.(l), expressed now using the Hubble function 
H as: 



8ttG 



c 



2 



tf 2 ~p' + y A (34) 

Where A is considered as a function that could change with z. The most 
general variation could be, from Eq.(34), expressed as: 

A = $i? 2 + *// + A (35) 

Where $, W and Ao are considered constants. It is easy to show that the 
general expression (35) is equivalent to: 

A = aip' + a (36) 

Where ao and a\ are new constants that are function of the previous ones. 
Taking into account all what has been stated, the dark energy density can be 
expressed as: 
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PA = ^A = /V+A> (37) 

Where (3$ and (3\ are new constants to be fixed. Considering the present 
values for the dark energy and matter (baryonic and dark) density parameters 
it can be deduced that ao = and ot\ — 5< ^f , so Eq.(36) can be written as (in 
agreement with [13]): 



And Eq.(37) will be: 



. 56ttG , 

A = ^ (38) 



PA - /V = \p' (39) 
Taking p = p' + p\ and with (39) it can be deduced: 

"-ITS' -55' (40) 

Eq.(40), and using (19), gives the same z dependence for the dark energy 
density parameter as the one it was deduced in case 4.2.2. The constant /3i 
could be let as a free parameter but we have taken the value deduced from (38) 
and (39) to use later. 

In this case the Eq.(7) will be the sum of the next pressures: 

p = wp=p m +p r +p A +Pi (41) 
Considering (40) and that the matter's pressure is zero, it can be deduced: 

(w + 0.7)p = Pr + Pl (42) 
Using (18), (19) and (26) in (42) we can deduce the interaction pressure: 

= i^Hl (z + f3\ 3 (a + z) (l {l + z) 0-1 \ 
ll 8irG \ 13 ) a \3 (a + z) f3 + z J 

3c 2 Hg ( 0.0001 4 



8ttG 



(1 + zT (43) 



Fig. 11 shows the interaction pressure which only reaches negative values on 
the future, a different behavior that showed in Fig. 8. 
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Figure 11: Interaction pressure for z running from -1 to 2. 

Finally Fig. 12 shows the dark energy density (40), which always takes posi- 
tive values. 
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Figure 12: Dark energy density for z running from -1 to 5. 

5 Conclusions. 

In this paper it has been shown how using a model for the universe, where 
the interaction between their components is admitted, it is possible to fit an 
effective EOS with a result similar to the models that not consider interaction 
(so assuming that the total pressure is only the sum of the pressures of each 
component). The studied EOS for the fluid of the universe varies with time and 
has three parameters that have been calculated using the SNelA experimental 
data. The final effective EOS for the universe fits the EOS for radiation for the 
early universe and at present time allows an accelerated era. 

This general model provides several signatures that could be used to see its 
suitability. Taking Hq = lOOh for the Hubble constant, the best value found 
for h has been 0.697. The predicted red shift, for the start of the actual phase 
of acceleration, is 2 £ [1.6, 1.7]. From Fig. 7 it is noticed that the model also 
predicts a maximum for - for a red shift z £ [0.31,0.33]. 
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From the general model it is possible consider different sub-models, but each 
one will always fulfill the above properties. In the paper it has been considered 
three sub- models that are physically relevant: two models with interaction, one 
for constant dark energy density and pressure and the other for variable dark 
energy density and pressure. The third model studied was without interaction 
and variable dark energy pressure, including two sub-cases. Some interesting 
conclusions has been deduced from the above models: the lack of necessity to 
consider phantom energy in order to fit a suitable model for the SNela ex- 
perimental data; On the other hand, the necessity of interaction between the 
different components of the universe to get positive dark energy density. Finally 
it has also been shown that there could be some kind of interaction in the dark 
sector (between dark energy and dark matter) and also between the dark sector 
and the baryonic matter. 
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